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I outline the power counting scheme recently introduced by M. Savage, M. Wise 
and myself for the effective field theory treatment of NN scattering. It is particu- 
larly useful for describing systems with a large scattering length, and differs from 
Weinberg's power counting. A renormalization group analysis plays a big role in 
determining the order of a given operator. Pions are ignored in this discussion; 
how to incorporate them is discussed in M. Savage's talkEJ. 



1 Introduction 

The purpose of this talk is to outhne an effective field theory expansion for 
NN scattering recently developed in Ref. □ with my coUaborators When 
pions are included (Martin Savage's talk in this volume), this expansion dif- 
fers from the original "jDnVO'jfi^ fj'^iY^n '^iV'^ Weinberg u and implemented 
or discussed since Refs. Q'Q'Elffl'E^rEjOllj'Lj. The advantage we claim for our 
scheme is that it allows analytic calculations of two-nucleon processes in a con- 
trolled expansion. As a result, the calculations of low energy NN processes 
are on the same footing as meson interactions in chiral perturbation theory 
(xPT) , or meson-baryon interactions in heavy-baryon chiral perturbation the- 
ory (HBxPT) . An obvious advantage of this approach is that it is much easier 
to extract physics from an analytical formula than from a computer simula- 
tion; for example, when one computes a Green function in terms of external 
momenta from Feynman diagrams, on has the whole analytic structure over 
a wide kinematic range and can analyze inelastic processes by examining the 
cuts, etc. Furthermore, one can better dissect the result and understand the 
relative importaiice from the various contributions. Finally, as we demonstrate 
explicitly in Ref.u, by being able to perform analytic calculations one can retire 
the old spectra of "off-shell ambiguity" , showing it to be hardly worth the fear 
and loathing that it usually inspires. No S'-matrix element depends on off-shell 
quantities! 



"Talk presented at the Joint Caltech/INT Workshop on Nuclear Physics with Effective Field 
Theory, February 1998. Report No. DOE/ER/40561-6-INT98 

''These notes contain somewhat more than my talk did, not out of revisionist tendencies, but 
because discussions at the workshop spurred us to compute the electromagnetic form factors 
of the deuteronQ, and in the process I have expanded my understanding of the subject. 
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One might think that in choosing to do an effective field theory calcula- 
tion of NN interactions analytically rather than numerically, one is sacrificing 
accuracy for elegance. In fact the power counting suggests this is not the 
case... Weinberg's proposal to expand the potential, and then resum its effects 
to all orders via the Lippmann-Schwinger equation to get the S'-matrix is not 
more accurate than simply expanding the S'-matrix directly. Furthermore, it 
does not even constitute an expansion unless one chooses the cutoff in the 
theory with some care. Therefore I see no virtue and some disadvantages with 
performing an effective field theory expansion of the potential, and then solving 
the Lippmann-Schwinger equation numerically. I believe that nuclear physics 
could well profit from abandoning such traditional tools as potentials and wave 
functions. 

Much emphasis has been placed on regulating and renormalizing the effec- 
tive theory of NN interactions. In fact, I view this workshop as the opportunity 
to put this issue behind us. The only reason to focus on these issues is to de- 
velop a consistent power counting scheme; once such a scheme is developed 
and followed consistently, results will not be scheme dependent, or else one 
is working with a model of NN interactions, as opposed to an effective field 
theory. So, while I emphasize in this talk the renormalization scheme used, 
the purpose is to develop the tools so that we can do physics, and not out of 
a love for formalism. 

2 Effective field theory for nonrelativistic scattering: a toy example 

Consider a toy model of heavy spinless "nucleons" N interacting via a Yukawa 
interaction characterized by a scale A. For low energy scattering we can con- 
struct the effective field theory describing scattering at momenta p <C A, con- 
sisting entirely of contact interactions in a derivative expansion. Since these 
local operators are singular, the formulation of the low energy theory neces- 
sarily introduces divergences that can be dealt with by conventional regular- 
ization and renormalization procedures, so that the final result is independent 
of a momentum cutoff. I will show how to organize the Feynman graphs in the 
effective theory in a consistent power counting scheme so that the scattering 
amplitude can be expanded in powers of p/A. Since the sizes of all the cou- 
pling constants in the effective theory depend on the subtraction scheme used 
to render diagrams finite, the development of the power counting scheme is 
intimately related to the renormalization procedure used. I will explain why 
the PL'S' subtraction scheme introduced in Ref.13 is particularly well suited for 
this problem. 

It should be no surprise that the effective field theory expansion for the 
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toy system is simply related to the conventional effective range expansion, and 
so the machinery of quantum field theory may appear to be heavy handed 
and superfluous. Nevertheless, the field theoretic language that I develop in 
this section is readily extended to the realistic problem of interest: nucleons 
interacting via both short-range interactions and long range pion exchange 
(see M. Savage's talk. Rem). In the realistic problem, effective field theory is 
not equivalent to an effective range expansion, and is the only framework that 
can consistently incorporate chiral symmetry and relativistic effects without 
resorting to phenomenological models. 

Assume that the spinless bosons N are nonrelativistic with mass AI, carry 
a conserved charge ( "baryon number" ) , and interact via the exchange of a 
meson cj) with mass A and coupling g. At tree level, meson exchange gives rise 
to the Yukawa interaction 

nO = -f^, (1) 
47r r 

and the Schrodinger equation for this system may be written as 

(2) 



X 



A2 



fEEAf, V^^, p'^ME. (3) 

Note that p is the magnitude of the 3-momentum carried by each N particle 
in the center of mass frame. Evidently there are two options for a perturbative 
solution for the S'-matrix for this system. The first is an expansion in powers of 
77, the familiar Born expansion. An alternative is to expand in powers of p/A, 
which is the expansion parameter used in effective field theory. An important 
feature of the low energy expansion is that it can provide accurate results in 
terms of a few phenomenological parameters even for nonperturbative rj. I will 
assume throughout that r/ ~ 1, since this is the regime we are interested in (a 
strongly coupled system without a plethora of bound states). 

The quantity that is natural to calculate in a field theory is the sum of 
Feynman graphs, which gives the amplitude iA, related to the S'-matrix by 

S^l + ^^A. (4) 
For S-wave scattering, A is related to the phase shift S by 

M p cot 5 ~ ip 
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From quantum mechanics it is well known that it is not A, but rather the 
quantity p cot S, which has a nice momentum expansion for p ^ A (the effective 
range expansion): 



p cot S = — 




where a is the scattering length, and rg is the effective range. So long as 77 ^ 1 
the coefficients r„ are generally 0(1/A) for all n. However, a can take on any 
value, diverging as rj approaches one of the critical couplings rjk for which there 
is a boundstate at threshold. (The lowest critical coupling is found numerically 
to be rji = 1.7.) Therefore the radius of convergence of a momentum expansion 
of depends on the size of the scattering length a. First I consider the situation 
where the scattering length is of natural size \a\ ~ 1/A, and then I discuss the 
case \a\ ^ 1/A, which is relevant for realistic TV TV scattering. 



2. 1 The momentum expansion for a scattering length of natural size 

In the regime \a\ 1/A and |r„| ~ 1/A, the amplitude A has a simple mo 
mentum expansion in terms of the low energy scattering data, 

Ana 



A 



M 



[1 - lap + (aro/2 - a^)p^ + O(pVA')] , 



(7) 



which converges up to momenta p ^ A. It is this expansion that we wish to 
reproduce in an effective field theory. 

The effective field theory of N particles interacting through contact inter- 
actions has the following Lagrangian: 



Ceff = {idt + V^/2M) N 



4-D 



Co(TVtiV)2 + ^ {NN)\N\/^N) + h.c 



,(8) 



where 



2V • V + V 



(9) 



The sum of Feynman diagrams computed in this theory gives us the amplitude 
A. As I will be using dimensional regularization for the loop integrals in 
this theory, the spacetime dimension is given by D ^. The ellipsis indicates 

'^Dimensional regularization is the preferred regularization scheme as it preserves gauge sym- 
metry and chiral symmetry, as well as Lorentz invariance (or Galilean invariance, for nonrel- 
ativistic systems). The advantage of the latter is that it makes the Feynman integrals easier 
to perform, as one can shift the integration variable. 
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Figure 1: The buhhle chain arising from local operators. The vertex is given by the tree level 
amplitude, eq. old). 



higher derivative operators, and is an arbitrary mass scale introduced 

to aUow the couphngs C2„ multiplying operators containing V^" to have the 
same dimension for any D. I focus on the s-wave channel (generalization to 
higher partial waves is straightforward), and assume that M is very large so 
that relativistic effects can be ignored. The form of the C2 operator is fixed by 
Galilean invariance, which implies that when all particle momenta are boosted 
p — > p + Mv, the Lagrangian must remain invariant. There exists another two 
derivative operator for p-wave scattering which I will not be discussing. 

In general, the tree level s partial wave amplitude in the center of mass 
frame arising from C^f f is 

00 

zA[:TJ = E ^^2„(a^)/" , (10) 

where the coefficients C2n (/^) are the couplings in the Lagrangian of operators 
with 2n gradients contributing to s-wave scattering. One may always trade 
time derivatives for spatial gradients, using the equations of motion when com- 
puting 5-matrix elements, and so we ignore such operators (see appendix B in 
RefQ for an explicit example).. 

The loop integrals one encounters in diagrams shown in Fig. ^ are of the 
form 

J {2n)D (^j^i^ + _ ill + -q^-^+^e) 
" ' J (27r)(^-i)^ \E-ceiM + ie) 

= -MiMEn-ME - ^er-^y-T (^) i^^^ . (11) 
In order to define the theory, one must specify a subtraction scheme; dif- 
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ferent subtraction schemes amount to a reshuffling between contributions from 
the vertices and contributions from the the UV part of the loop integration. 
How does one choose a subtraction scheme that is useful? I am considering 
the case \a\, |r„| ~ 1/A, and wish to reproduce the expansion of the amphtude 
eq. (|^). In order to do this via Feynman diagrams, it is convenient if any Feyn- 
man graph with a particular set of operators at the vertices only contributes to 
the expansion of the amplitude at a particular order. Since the the expansion 
eq. is a strict Taylor expansion in p, it is it is therefore very convenient if 
each Feynman graph gives one a simple monomial in p. Obviously, this won't 
be true in a random subtraction scheme. A subtraction scheme that fulfills 
this criterion is the minimal subtraction scheme (MS) which amounts to sub- 
tracting any 1/(1? — 4) pole before taking the D ^ A limit. As the integral 
eq. ( pj] ) doesn't exhibit any such poles, the result in AIS is simply 

tMS _ /n,rj:^\n f M\ , — j- _ ,/Af\ 2,H 



/,f ^ = {MEr (^^ j V-ME-te = j p^^+i . (12) 

Note the nice feature of this scheme that the factors of q inside the loop get 
converted to factors of p, the external momentum. Similarly, a factor of the 
equations of motion, idt + V^/2M, acting on one of the internal legs at the 
vertex, causes the loop integral to vanish. Therefore one can use the on-shell, 
tree level amplitude eq. (|l^) as the internal vertices in loop diagrams; summing 
the bubble diagrams in the center of mass frame gives 



EC2nP^ 

1 + liMp/Air) J: C72„p2 



-4 = ...tr.tr^.. (13) 



Since there are no poles at _D = 4 in the MS scheme, the coefficients C2n 
are independent of the subtraction point fi. The power counting in the MS 
scheme is particularly simple, as promised: 

1. Each propagator counts as 1/p^; 

2. Each loop integration J dL^q counts as p^ (since go ~ q^/2Af); 

3. Each vertex C2„V^" contributes p^". 

The amplitude may be expanded in powers of p as 

00 

A = J2An, An^Oipn (14) 

where the An each arise from graphs with L < n loops and can be equated 
to the low energy scattering data eq. (m) in order to fit the C2n couplings. In 
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particular, Aq arises from the tree graph with Co at the vertex; Ai is given by 
the 1-loop diagram with two Cq vertices; A2 is gets contributions from both 
the 2-loop diagram with three Cq vertices, as well as the tree diagram with 
one C2 vertex, and so forth. Thus the first three terms are 

A = -Co, A^^^C'Q^, A,^C'o[^)' -C2P' ■ (15) 

Comparing eqs. (0, |l5|) I find for the first two couplings of the effective theory 

C„ = ^, C.=Co^. (16) 
In general, when the scattering length has natural size, 

^'"^]^A^- ^^^^ 

Note that the effective field theory calculation in this scheme is completely 
perturbative even though ~ 1 and there may be a boundstate well below 
threshold. The point is, that when there are no poles in A in the region |p| <^ A, 
the amplitude is amenable to a Taylor expansion in p/A in that region; with a 
suitable subtraction scheme (MS), this Taylor expansion can correspond to a 
perturbative sum of Feynman graphs. 

2.2 The momentum expansion for large scattering length 

Now consider the case \a\ 3> 1/A, |r„| 1/A, which is of relevance to realistic 
NN scattering. For a nonperturbative interaction (77 ~ 1) with a boundstate 
near threshold, the expansion of A in powers of p is of little practical value, 
as it breaks down for momenta p > l/\a\, far below A. In the above effective 
theory, this occurs because the couplings C2n sue anomalously large, C2n ~ 
47ra"+^/A/A". However, the problem is not with the effective field theory 
method, but rather with the subtraction scheme chosen. 

Instead of reproducing the expansion of the amplitude shown in eq. (0), 
one needs to expand in powers of p/A while retaining ap to all orders: 



47r 1 

A 



M (l/a + ip) 



rV2_ 2 , {ro/^r 4 , (^i/2A^) 4 
{1/a + ipf ^(l/a + ip)2^ ^{1/a + ipf 



(18) 



Note that for p > l/\a\ the terms in this expansion scale as {p . . .}. 

Therefore, the expansion in the effective theory should take the form 



^ - ^ , An- 0{p'') (19) 
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beginning at n = —1 instead of n 
with eq. ([1^), we see that 

A-i = 

Aq = 

and so forth. Again, the task is to compute the An in the effective theory, 
and equate to the appropriate expression above, thereby fixing the C2n coeffi- 
cients. As before, the goal is actuaUy more ambitious: each particular graph 
contributing to An should be 0{p"), so that the power counting is transparent. 

As any single diagram in the effective theory is proportional to positive 
powers of p, computing the leading term A^i must involve summing an infinite 
set of diagrams. It is easy to see that the leading tei-m A-i can be reproduced 
by the sum of bubble diagrams with Co vertices 0, which yields in the AIS 
scheme 

Comparing this with cq. ( [20| ) gives Co — Ana/M, as in the previous section. 
However, there is no expansion parameter that justifies this summation: each 
individual graph in the bubble sum goes as Co(CoMp)^ ^ (47ra/M)(iap)^, 
where L is the number of loops. Therefore each graph in the bubble sum is 
bigger than the preceding one, for \ap\ > 1, while they sum up to something 
small. 

This is an unpleasant situation for an effective field theory; it is impor- 
tant to have an expansion parameter so that one can identify the order of any 
particular graph, and sum the graphs consistently. Without such an expan- 
sion parameter, one cannot determine the size of omitted contributions, and 
one can end up retaining certain graphs while dropping operators needed to 
renormalize those graphs. This results in a model-dependent description of the 
short distance physics, as opposed to a proper effective field theory calculation. 

Since the sizes of the contact interactions depend on the renormalization 
scheme one uses, the task becomes one of identifying the appropriate subtrac- 
tion scheme that makes the power counting simple and manifest. The MS 
scheme fails on this point; however this is not a problem with dimensional reg- 
ularization, but rather a problem with the minimal subtraction scheme itself. 
The momentum space subtraction at threshold used in Ref.H behaves similarly. 

Next, consider an alternative regularization and renormalization scheme, 
namely to using a momentum cutoff equal to A. Then for large a one finds Co ~ 



0, as in the expansion eq. (|Tj). Comparing 
47r 1 

M(l/a + ip)2 ' ^ 
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(47r/MA), and each additional loop contributes a factor of Co (A + ip)M /Att ~ 
(1 + ip/A). The problem with this scheme is that for A ^ p the term ip/A 
from the loop is small relative to the 1, and ought to be ignorable; however, 
neglecting it would fail to reproduce the desired result eq. (^^. This scheme 
suffers from significant cancellations between terms, and so once again the 
power counting is not manifest. 

Evidently, since A-i scales as l/p, the desired expansion would have each 
individual graph contributing to A-i scale as l/p. As the tree level contri- 
bution is Co, I must therefore have Co be of size cx l/p, and each additional 
loop must be 0(1). This can be achieved by using dimensional regularization 
and the PDS (power divergence subtraction) scheme introduced in Ref . □. The 
PDS scheme involves subtracting from the dimensionally regulated loop inte- 
grals not only the 1/(1? — 4) poles corresponding to log divergences, as in MS, 
but also poles in lower dimension which correspond to power law divergences 
at D = 4. The integral /„ in cq. ( [Til ) has a pole in. D — 3 dimensions which 
can be removed by adding to /„ the counterterm 

so that the subtracted integral in D = A dimensions is 

/^''^ = In + SIn = - {MET (^) (a* + ■ (23) 



In this subtraction scheme 



47r 



47r 



+ IJ. + ip 



(24) 



By performing a Taylor expansion of the denominator of the above expression, 
and comparing with eq. (p^, one finds that for fj. 3> l/|a|, the couplings C2n(/^) 
scale as 

CM - . (25) 

Eq. p5| implies that fJ, ^ p, C2„(^) ~ l/p""*"^. A factor of V^" at a vertex scales 
as p , while each loop contributes a factor of p. The power counting rules for 
the case of large scattering length are therefore: 

1. Each propagator counts as 1/p^; 

2. Each loop integration J d^q counts as p^; 
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3. Each vertex C2nV^" contributes p" 



We see that this scheme avoids the problems encountered with the choices of 
the MS {fi = 0) or momentum cuttoff (/i ~ A) schemes. First of aU, a tree 
level diagram with a Cq vertex is 0{p~'^), while each loop with a Cq vertex 
contributes Cq{p)M{h + ip)/47r '-^ 1. Therefore each term in the bubble sum 
contributing to A^i is of order p^^, unlike the case for = 0. Secondly, since 
/X ~ p, it makes sense keeping both the jjL and the ip in eq. (23) as they are 
of similar size, unlike what we found in the fj, = A case. The PDS scheme 
retains the nice feature of MS that powers of q inside the loop integration are 
effectively replaced by powers of the external momentum p |^. 

Starting from the above counting rules one finds that the leading order 
contribution to the scattering amplitude A-i scales as p~^ and consists of 
the sum of bubble diagrams with Co vertices; contributions to the amplitude 
scaling as higher powers of p come from perturbative insertions of derivative 
interactions, dressed to all orders by Cq. The first three terms in the expansion 
are 



A-i = 
Ao = 

Ai = 



-Co 



[1 



-C2P^ 



ip)] 



( (C2p2)2M(M + ip)/47r 



[l + ^{t. + ^p) 



[1 



CqM 
4tt 



(P + ip)] ' 



(26) 



where the first two correspond to the Feynman diagrams in Fig. Q. The third 
term, Ai, comes from graphs with either one insertion of C4V'* or two insertions 
of C2V^, dressed to all orders by the Co interaction. 

Comparing eq. ( |2^ ) with the expansion of the amplitude eq. (^8|), the 
couplings C2n are related to the low energy scattering data a, r„: 



C^o(/i) = 

C2i^i) = 
Cii^l) = 



An 
M 
An 
M 

An 
M 



1 



—pi + 1/a 



1 



IE 

\/a) 2 

3 



+ 1/a 



1 



1 ri 

2 A2 



1/a) 



(27) 



^An alternative subtraction scheme with similar powcrre«unting is to perform a momentum 



subtraction at p 



as recently suggested in Ref. m. 
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Figure 2: Leading and subleading contributions arising from local operators. 



Note that assuming r„ ~ 1/A, these expressions are consistent with the scahng 
law in eq. (|25|). 



2.3 The renormalization group 

This power counting described in the previous section relies entirely on the 
running of C2n(/i) as a function of /i given in eq. (|2^). This was derived by 
summing up all of the diagrams in Fig. ^ explicitly, and then comparing with 
the form of a general amplitude, eq. (pf). When pions are included in real NN 
interactions, the diagrams in Fig. ^ cannot be explicitly summed. However, the 
power counting can be established perturbatively by examining the /3-functions 
and the renormalization group running of the C2n(M) couplings. The depen- 
dence of C2n (m) on /i is determined by the requirement that the amplitude be 
independent of the arbitrary parameter /i. The physical parameters a, rn enter 
as boundary conditions on the RG equations. 
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ft 



Fi gure 3: Graphs contributing to the (^-functions for C2n 



The /^-function /?2n for the couphng C271 is defined by 

dC2n 



d/i 



(28) 



and all of the /3-functions can be computed by requiring that any physical 
quantity (e.g. the scattering amplitude) be independent of /i. In the PDS 
scheme, the dependence of the C'2n coefficients enters either logarithmically 
or linearly, associated with simple 1/(0 — 4) or 1/{D — 3) poles respectively. 
The functions f32n follow straightforwardly from /i^(l/^) = 0, using the 
expression for A in eq. (M). This gives 



3 -^^/^V 
- ^ 2. 



C2mC'2(n-m) 



(29) 



However, I did not need the full, explicit amplitude A to compute (32n, as 
the exact /3-functions can be computed from the one-loop diagrams shown in 
Fig. ^. That is because only the one-loop diagrams contribute to simple poles. 
Without pions, the only poles one encounters are all of the form 1/(_D — 3). 
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We examine the RG equations for the first two couphngs, Co and C2, in 
order to expHcitly show how one recovers the resuhs in eq. ( p7| ) from solving the 
renormahzation group equations. From eq. ( p9| ) — or equivalently, the diagrams 
in Fig. ^ — it foUows that 

/32 = 2^CoC2 . (30) 
47r 

Integrating these equations relates the C2n coefficients at two diS'erent renor- 
mahzation scales /i and /xq. Comparing the theory with A and its derivatives 
at /i = p = determines the initial values C2n (0) as in eq. ( |l6| ) . The solutions 
to eq. ( ^0|) are 



l + Co(Aio)M(Mo-A*)/47r ' 



2 



The boundary conditions are supplied by equating the computed A to scatter- 
ing data at = p = 1^ which yields 

Co(0) = 4™/M , C2(0) = Co(0)aro/2 , (32) 

and so forth. With these boundary conditions and the solutions eq. (^) we 
arrive at the result derived previously for Co(a*) ^'iid C2(/x) in eq. (p7|): 



It is instructive to solve the complete, coupled RG equation 

IJ'-T-C2n — ^ C2mG2{n~m) (34) 

m=0 



'^This is not the only choice. For example, when discussing the deuteron, it is convenient to 
fix the location of the pole to the true deuteron binding energy. The differences arising from 
this alternate fitting procedure are small higher order effects. 
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for the leading small fi behavior of each of the coefficients C2n ■ The solution, 
for n > is 

M(-^ + l/a) \-fi+ 1/a J 

First note that the scaling property in eq. ( p5| ) is realized: C2„(/i) oc 
for \l/a\ <C n <C A. What is curious is that this leading behavior does not 
entail a new integration constant for each n, but only depends on the two 
parameters a and ro encountered when solving for Co(/i) and C2(/i); this is 
due to a quasi-fixed point behavior of the RG equations — the C2n couplings 
are being driven primarily by lower dimensional interactions. One can see this 
explicitly in our formula eq. (|2^) for C4, where the leading 0{p,~^) part of C4 
depends only on ro, while the subleading 0{ii~^) part is proportional to ri. 

This behavior allows us to establish a connection between the prese]at|Work, 
and the method of introducing an s-channel dibaryon discussed in Ref. c3. The 
leading behavior of all of the C2n coefficients is determined by the effective 
range ro . If one resums this leading behavior at the NN vertex one finds (for 
Ai-p> l/|a|) 



A ^ 

(Stt/MVo) 



E~ {~^l+\/a)/Mro 



(36) 



This looks like an s-channel propagator for a particle at rest of mass [2M + 
{—fj, + l/a)/Mrol— and in fact, for /j, — 0, corresponds exactly to the dibaryon 
proposed in Ref. E3 to reproduce the scattering due to a short range potential. 
We see that using the dibaryon is as good as (but no better than) carrying 
out the effective field theory calculation to O(p^). The subleading corrections 
can be accounted for by including the subleading part of the Ci{^)p^ vertex 
proportional to ri, and which occurs at 0{p^). Tljis dibaryon was recently 
used with great success in the three-body problemll3. 



3 Expanding the potential or the scattering amplitude? 

Throughout this talk I have discussed an expansion of the scattering ampli- 
tude, while most previous work in the subject of effec.tiYP field t heor y f or NN 
scattering has emphasized expansion of the potentia0^Bffl^BElil3'M'E3E3o, fol- 
lowed by a solution of the Lippmann-Schwinger equation with this approximate 
potential. Since every observable is an S-matrix element, it would seem that 
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the two methods, if carried out to the same order, ought to give answers that 
disagree only by higher order effects. For example, if working to 0(jP), the di- 
rect expansion of ~ A-i + Aq is linear in the C2V^ operator. Alternatively, 
if the potential is expanded to linear order in C2V^, the Lippmann-Schwinger 
equation yields an amplitude which reproduces the term linear in C2 , but also 
includes terms higher order in C2. These terms nonlinear in C2 would appear 
to be higher order effects and negligible. 

This is not in general true, however. Unlike the A expansion I have out- 
lined which is explicitly scheme independent, the Lippmann-Schwinger ap- 
proach is not. That is because the time ordered product of several insertions 
of the potential induces new divergences that require counterterms not in- 
cluded in the expansion. For example, when the potential is taken to include 
only the Cq and C2V^ operators, the one-loop contribution to the Lippmann- 
Schwinger equation has a divergence that requires the operator to absorb 
the divergence — an operator not included in the expansion! As a result, the 
Lippmann-Schwinger approach is necessarily scheme-dependent. 

Is this bad? Not fatal, but undesirable. It means that the size of neglected 
effects depends on the renormalization scheme and cutoff or renormalization 
scale. That is why in calculations of this sort, the renormalization scale or 
cutoff becomes an extra parameter that has to be chosen to minimize the 
errors. 

Understanding the scheme dependence of the Lippmann-Schwinger result 
helps explain an apparent paradox: the amplitude A calculated in this talk 
in the PDS scheme is /i independent at each order in the expansion. In 
particular, each term An is unchanged if I take — > 0, which is the MS 
scheme. However, there have been numerous discussions about how MS is sick. 
In fact, MS is not sick except for the fact that at any given order there tend 
to be large cancellations among the graphs included at that order. However, 
in the Lippmann-Schwinger approach, parts of higher order terms are kept, 
which destroys the cancellations; the result is fj, dependent, and the expansion 
becomes very bad for fi <^ p. For example, at 0{p^), Ai given in eq. (26) 
involves both a C| term and a C4 term; for ji p these two terms are both 
large and cancel against each other. However the 0{p^) Lippmann-Schwinger 
calculation, in which the potential includes the C2 interaction, but not the C4 
term, gives a scattering amplitude that includes the C| piece of Ai but not the 
C4 piece which is needed to largely cancel the C| contribution and render its 
effects small. If one solves the Lippmann-Schwinger equation and chooses the 
PDS value ^ ^ p, then the C| and C4 terms in Ai are of the same size and (by 
construction) and small, and the error in the Lippmann-Schwinger apiplitude 
really is higher order. This has been independently discussed in Ref. 
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The correct statement about the MS subtraction scheme is that it can 
be misleading, in that one cannot look at an individual Feynman graph at 
contributing to An and expect it to contrihuta^ 0{p'^). The problems with 
the MS scheme reported in the literature LraHo are the result of combining 
the MS scheme with the scheme dependent Lippmann-Schwinger approach. I 
favor avoiding the Lippmann-Schwinger approach altogether — why would one 
want a calculational scheme where the size of the errors at a given order in the 
expansion depend on the renormalization scheme? In effective field theory, with 
its necessarily singular interactions that require renormalization, the concept 
of a classical potential does not seem very useful. 



4 Conclusions and a challenge 

I hope I have convinced you that there exists a rather simple and well defined 
approach to calculating low energy processes involving two nucleons. There 
are quite a few interactions of interest to calculate, such as NN NNj, 
NN dj, pp — > de^v, nd —>■ nd, as well as isospin violation and parity 
violation in NN scattering. And for every calculation, there is always one 
higher order that can be computed! 

However, it would be disappointing if the applicability oLthese techniques 
were limited to two nucleon processes. Evidence to date clu suggests that 
the expansion works reasonably well up to momenta comparable to the Fermi 
momentum in nuclear matter. What are the prospects of using effective field 
theory to discuss the structure of nuclei, or the equation of state of nuclear 
matter? 

At present the road block is the three-body interaction. It is tempting to 
dismiss the three-body interaction as negligible, as often claimed. However, 
the two-body interactions renormalize the three-body force, which means that 
in fact the strength of the three-body interaction is scheme dependent, and 
saying it is aniall in general is nonsensical. In fact, because of the EfimovEj 
and ThomasEj effects, we know that we cannot keep zero-range two-body in- 
teractions while neglecting three-body interactions. It is clear, due to Pauli 
statistics, that at most there can be four-body contact interactions (without 
derivatives). However, until the properties of the three- and four- body inter- 
actions are understood properly, it is unclear how to proceed to a discussion 
of nuclear matter. 1 view understanding this issue of few-body interactions to 
be the fundamental challenge in this field. 
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